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Roles of antinucleon degrees of freedom in the relativistic random phase approxima- 
tion(RPA) are investigated. The energy-weighted sum of the RPA transition strengths 
is expressed in terms of the double commutator between the excitation operator and 
the Hamiltonian, as in nonrelativistic models. The commutator, however, should not 
be calculated with a usual way in the local field theory, because, otherwise, the sum 
vanishes. The sum value obtained correctly from the commutator is infinite, owing to 
the Dirac sea. Most of the previous calculations takes into account only a part of the 
nucleon-antinucleon states, in order to avoid the divergence problems. As a result, 
RPA states with negative excitation energy appear, which make the sum value vanish. 
Moreover, disregarding the divergence changes the sign of nuclear interactions in the 
RPA equation which describes the coupling of the nucleon particle-hole states with 
the nucleon-antinucleon states. Indeed, excitation energies of the spurious state and 
giant monopole states in the no-sea approximation are dominated by those unphysi¬ 
cal changes. The baryon current conservation can be described without touching the 
divergence problems. A schematic model with separable interactions is presented, 
which makes the structure of the relativistic RPA transparent. 


1 Introduction 

It has been shown by many authors [UElEJIl] that relativistic models, assuming nuclei to be 
composed of Dirac particles and various mesons, work well phenomenologically to reproduce 
nuclear static and dynamic properties. In most of calculations, however, antinucleon(N) 
degrees of freedom are not fully taken into account, in spite of the fact that they are one of 
the main differences between relativistic and non-relativistic models. The reason why the 
full space is not used is because there exist divergence problems [H [5l [6] which are not yet 
handled with a proper way for finite nuclei. 

Although all the space was not included in the previous calculations, a part of N degrees 
of freedom was taken into account, aiming to keep some fundamental principles and to 
reproduce gross properties of nuclei. 

For example, in the random phase approximation (RPA), the baryon current conserva¬ 
tion required some excitations of antinucleons O [H [7]. It was necessary for description of 
the center of mass motion to have the Landau-Migdal parameter Fi taking into account a 
part of N(nucleon)-N states[71[8]. In the same way, the spurious states in RPAf^fTO] was de¬ 
scribed with those N-N excitations in addition to nucleon particle-hole states. Furthermore, 
an abnormal enhancement of the isoscalar magnetic moment in the Hartree approximation 
demanded corrections from the N-N excitations through Fi m- 

For the above reasons, the two ways to include the N space, avoiding the divergence 
problem, were proposed. The one is the no free term approximation(NFA) which simply 
neglects the divergent part of the RPA response functions. The remaining part is composed 
of transitions of antinucleons to the Fermi sea which are in fact Pauli-blocked [Sj [6]. The 


1 


other is the no-sea approximation (NS A) to assume that all the N states are empty. There is 
no divergence problem in this case also, but transitions of nucleons in the Fermi sea to the 
N states are permitted with negative excitation energies|10j. Even though both methods 
look unreasonable, they have been widely used, since experimental data are well reproduced 
phenomenologically [3l [T^ 1131 [H] . In fact, as will be seen later explicitly, NFA and NSA are 
equivalent to each other in RPA|10j. 

The purpose of the present paper is to investigate the structure of the relativistic RPA 
in detail, and in particular, a role of N degrees of freedom there. The relativistic RPA can 
be developed almost in the same way as non-relativistic ones. It will be shown, however, 
that a careful treatment is required for N degrees of freedom. They cause the divergence 
problems, but cannot be simply ignored as in NFA and NSA. 

It will be shown that the energy-weighted sum of the RPA transition strengths is ex¬ 
pressed formally in terms of the double commutator between the transition operator and 
the Hamiltonian. This RPA sum rule is the same as in nonrelativistic case. The commu¬ 
tator, however, should not be calculated employing a usual rule in the local field theory. 
Otherwise, the sum vanishes, in spite of the fact that it should be positive definite. The 
correct calculation of the commutator gives the sum value to be infinite, because of N-N 
excitations. 

In NFA and NSA, RPA states with negative excitation-energy appear, in addition to 
those with positive energy. This is due to the fact that the full N-N excitations are not take 
into account. As a result, the energy-weighted sum value of the excitation strengths vanishes. 
Moreover, disregard of the divergence in NSA and NFA gives rise to other unphysical results. 
The sign of the nuclear interactions is changed in RPA equation relevant to the remaining 
N-N states. Because of this fact, attractive (repulsive) forces work as repulsive (attractive) 
ones in the coupling of nucleon particle-hole states with N-N states considered in NFA and 
NSA. These effects are not negligible, but dominate the excitation energies of some low 
lying states. For example, the previous numerical calculations could reproduce a spurious 
stateOHU] and giant monopole states [121 [HI HI], invoking those unphysical effects hidden 
in the RPA equation. 

The above defects of NFA and NSA have not been investigated explicitly so far, since 
the structure of the relativistic RPA formulae may be rather complicated, compared with 
non-relativistic ones. In the next section, we will briefly review the relativistic RPA focusing 
on roles of N-N states in various approximations. In ^ the energy-weighted sum of the 
RPA excitation strengths is discussed, according to recent new insight for the long-standing 
problem on the relativistic sum rule|15]. In 21 21 will explore in detail why NFA 

and NSA seemed to describe well the continuity equation and the spurious state without 
the full N-N excitations. 

It is well-known that a schematic model with separable interactions helps us to under¬ 
stand the structure of nonrelativistic RPA|16]. In 21 3. similar model will be introduced 
for the relativistic RPA, in order to make clear the present discussions. The final section is 
devoted to a brief summary of the present paper. 

2 Relativistic RPA 

The relativistic RPA was formulated in various ways[Il El El El HO]. Taking notice of the 
dependence on N-N states, let us briefly review it. 

We begin with the Hartree polarization function for arbitrary 4x4 matrices, A and B 
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[S], 

inn{A{l),B{2)) = { lT(^(l)^V^(l)^(2)i?V'(2)) | ), (1) 

defined for the Hartree ground state | ). The baryon field ^(x) is written using the complete 
set of the eigenfunction, (pa, of Hartree Hamiltonian, 

= '^(Pa{x)e~"^°‘^aa, (2) 

a 

where denotes the eigenvalue of and plays a role of the annihilation or creation 
operator of N and N, satisfying { Oq , } = dap, etc.. Then, together with the closure 

property of (pa, 

‘fa{x)‘pi{x') = 6{x - x'), (3) 

a 

the baryon field satisfies the anti-commutation relation, 

{lpa{xo,x) , ll:l{xo,x') } = 6abS{x - x), ( 4 ) 

where a and b denote the Dirac matrix indices. Since the simple interactions have rather 
advantage for our purpose to investigate the detail of the relativistic RPA structure, we 
assume, throughout the present paper, the a-LO model[T], which provides us with the Hartree 
Hamiltonian as 


h{x) =-ia-V+ -foM+ Va{x), V}i{x) ='^^^{x). (5) 

In the above equation, we have employed the following abbreviation for the potential, 

E{x) = V„{x) + -i^^V^{x) (6) 

with 


y^{x)=gl I D^{x - y) { | V'(y)V’(y) I ), 

yti{x)=gl [d'^y D^{x - y) { \'tpiyht^'tp{y) \ ), 


(7) 

( 8 ) 


Qa and being the coupling constants, and and denoting the a- and cu-meson 
propagators, respectively [6], 


Da{x) = D„{ko,x) = j 

DUx)= f^e-^>‘°^°DUko,x), D^{ko,x) = - 


Ak-x 


d^k 


(27r)3 k^k>^ — ’ 


Akx 


d^k 


(27r)^ — rri^ + ie 


Here, tUo- and represent the masses of the a- and cj-mesons. 

The calculation of Eq.([T]) gives [6], 

nu{A{l),B{2)) = ^ Jdu;e-^‘^^^^-^^An{A{xi),B{x 2 ),uj), 

where we have defined 

77h(H(®i),H(®2),w) = Y^fi(^^)^^aixi)'^aix2)Bipp{x2) {nu)ap{<^) 

OL^ 


(9) 

( 10 ) 


( 11 ) 


( 12 ) 
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with 


{nR)al3{^) = 


Nap 


N, 


/3q 


Nap = (1 - 0c 


(13) 


oj — Eap + ie io — Eap — ie 

In this equation, Eap stands for the excitation energy, Ea — Ep, and Nap implies the 
transition from the occupied state to unoccupied state as. 


1, a : occupied state in | ), 

0 , a : unoccupied state in | ). 


(14) 


Defining the Fermi energy by Ep, it is convenient to employ the following notations for Nap, 

Pa = 0{Ea - Ep) , ha = 6{Ea) 0{Ep - Ea) , ha = 6{-Ea), (15) 

where pa, ha, and ha indicate the particle, hole and N states, respectively. In the full 
space calculation, the occupied states are expressed by {h + h), and unoccupied states by 
p, so that Nap is written as Nap = Pa{hp + hp). In NSA where the N states are empty, 
we have Nap = {pa + ha)hp. In writing Nap = Pahp — hahp + {pa + ha)hp for the full 
space, NFA neglects the last term which describes the vacuum polarization to produce the 
divergence[5l[6]. In the approximation neglecting all the N states(NoN), Nap is simply given 
hy Pahp. Thus, the difference between various approximations can be represented by Nap, 
reading it as 

Pahp + Pahp (Full), 

Pahp + hahp (NSA), 

Pahp-hahp (NFA), 

^ Pahp (NoN). 

-1 


NaP = < 


(16) 


We define the inverse of 77 h by iljj as 

{nYi^)ap,a"P"{npi)a"P",a'P' = IaP,a>P', 

a"P" 

where the following abbreviations are used 

(-^H)a/3,a'/3' = ^aa'^PP'{Nil)ap , Iap,a'p' = ^aa'^PP' (^q/3 + ' 

For = 1, we obtain 

(TTjj {ph))ap,a'P' — ^aa'^PP' {NaP Npa) {oJ Eap) ■ (f'^) 

Because of pa + /i^ + ha = 1, the above {Nap — Npa) is given for each approximation as 


NaP Nf^a — * 


(hp + 71 / 3 ) - {ha + ria) (Full), 
h0 - ha (NSA, NFA), 

Pahp-ppha (NoN). 


(18) 


Unlike iln, the inverse 77^^ of NSA is the same as for NFA. This fact makes NSA and NFA 
are equivalent to each other in RPA, as seen later. 

The RPA polarization function is written in terms of TIhIS], 

77(A(1),R(2))=7Ih(A(1),R(2)) 

+ gl Jd^xd^ynu{A{l),r\x))D,{x-y)n{r,{y),B{2)), (19) 
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where rj = —1 for a and, hence, = 1 for rj = — 1 , and 7 ^ for rj = jj,. As in TTh, we write 
Ea. (|19p in the form: 


n{A{l),B{2)) = ^ I dcae—(20 ) 

with 

n{A{xi),B{x2),uj) = 77h(A(®i),5(®2),w) 

+ ff^y d^xd^ynn{A{xi),r'^{x),uj)Dr^{u;,x- y)77(r^(y),5(®2),w). 

The above n{A{xi), B{x 2 ),uj) is described as 

n{A{xi),B{x2),U!) = Ffsixi)A(pa{xi) na0,a'y'{^^)Fa'ix2)B(f0'{x2), ( 21 ) 

0.^0.' ji' 

where we have defined 


nal3,a'l3'{x>) = 6aa'^PI3'{F[-R)aj3{^) + ^ {F[-R)ap{'^)Vap",pa"{^)F[a"p",a'P'{xi), 

a"P" 


( 22 ) 


using the notation 

Val3',0a'{x’) = 9^ J d^Xid^X2^a{xi)r'^ip^{xi)Drj{u},Xi - * 2 ) (a;2)-f'r,<^a'(* 2 ) • 

This is also written as 

Va/3',/3a'{^) = («,Vi 2 (a;) |/3,a'), (23) 

with 

Vi2(w) = 5^(7oT’')i(7or^)2-C>^(w,®i - X 2 ). 

In relativistic models, nuclear interactions contain cj-dependence coming from Eqs.Q and 
m- For later discussions of the relativistic RPA, however, we have to neglect this retar¬ 
dation effect, as in nonrelativistic RPA. Fortunately, their contributions to the interactions 
have been shown to be negligible in Ref. m, when ma,moj > w. Hence, from now on we 
will develop the RPA, assuming Vi 2 ( 0 ). 

Eg. (12211 is formally described as follows, 

n{u!) = 7Ih(w) + n}i{u)vn{uj) , (R'H)a/3,a'/3'(<^) = (5/3/3'(I?'H)a/ 3 (w). (24) 

Then, the inverse of F[(u) can be written as 7T~^ = — V with 

(7T {^')')alS,a'0' — (5qq.'( 5/3/3' (A^o/3 5V/3 q,) (cj Pa/3',/3a'' (25) 

Now let us define the eigenvector of the above eauation|18j. 

n-\un)c^'^'> =0, 

which is 

(5Va/3 - iV/3a) ^ Ra/3',/3a'= (A^a/3 " iV/3a) ' (26) 

q'/3' 

In writing 

(27) 
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the coupled equations are obtained, 


+ Vo.a'= ^nX^afl ^ 


a'p' 


+ Y. " UJnVY’■ 

a'p' 


Employing the abbreviations 


(28) 


-^q:/9,q^/3' — ^ olol '-|- NN Q^f^ 


Eq. (|28p provides us with the relativistic RPA equation of the form: 


E 

o'/3' 


-^o/3,a'/3' Bo^po^ipi 
^*ap,a'ji' ^al3,a'/3' 


^a'P' 

vl”'l 

^ a'P’ 


— CUn N(^P 


X 


(n) 

0.13 


-y 


(n) 


O^ 


(29) 


(30) 


When the N degrees of freedom in are neglected, the above equation rednces to the 
well-known non-relativistic RPA eauation |16l I19j . 

The relationship between the Fnll and NFA can be seen in Ea. (|27p . Comparing 
of NFA with the one of the Full case, the part of the N-N excitations in NFA has a mi¬ 
nus sign, as —hafi/s, because of neglecting the divergent part as mentioned before. This 
additional sign induces unphysical effects in NFA that attractive(repulsive) forces work as 
repulsive(attractive) ones in the /i^h^-dependent part of the RPA equation. Eg. ([281) . 

R is also seen from Ea. (|27l) that NS A is eqnivalent to NFA. In writing Ea. (|27l) explicitly 
for NS A and NFA as 

= {pahfj + nahf^)X^^^ + {ppK + 

Ca0^ = {pah^ - -k {pjiK - h0no)Yf^^, 

the replacements 

yNSA _ t^NFA wNSA _ _ y-NFA '^NSA _ -y^NFA y_NSA _ _ yNFA 

^ph ^ph ! ^nh ^ hn ’ ^ph ^ph > ^ nh ^hn 

lead to 

r«NSA _ r'NFA 

^ojB ^0/3 

Hence, NSA yields same unphysical change of the sign in the interactions, and the same 
eigenvalues, as in NFA. 

We add a few comments which we need for later discussions. First, the complex conjugate 
of Ea. (l26l) implies that cjj^* is also its solution with the eigenvalue to be —LUn- Second, 
the orthogonality and normalization of the eigenvectors are written as 

E - ^/3«) = ^nSnn' {K = ±1). (31) 

ajS 


The eigenvectors with = 1 describe the RPA excited states. Thus, the norm of the 
eigenvectors in NSA is also the same as in NFA. Third, the closure relation is given by 

E (Nap - iV/3a) ■ (32) 
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Finally, we note that Eas. pU|l and (IHTD provide us with 

y<”)t) M ( j , >1 = ( ^. ® ^ (33) 

In non-relativistic models, we usually obtain > 0 for A„ = 1[I6119]. In NS A and NFA, 
however, it is not necessary for A4 to be positive definite, so that there may be negative 
eigenvalues Un even for A„ = 1 . 

3 Energy-weighted sum of the excitation strengths 

One of the reason why RPA has extensively been used in non-relativistic nuclear models is 
because the sum of the energy-weighted strengths for exciting the RPA states is constrained 
in a simple way by the employed Hartree or Hartree-Fock basis. Once we have the double 
commutator of the one-body Hermitian operator F 

F = (34) 

with the Hamiltonian H, the sum is equal to its expectation value of the Hartree or Hartree- 
Fock ground state [T6l fl^ . 

j;a;„|(n|#|0)|2 = l( \[F,[H,F]]\ ). (35) 

n 

At first glance, however, it seems that the above non-relativistic rule is not extended to 
the relativistic RPA. If the operator F has the only coordinate-dependence, the double 
commutator with the Dirac equation vanishes, 

[F((r), [ho, F((r)]] =0, (36) 

since ho contains the only linear derivative, 

ho(x) = — foe V -|- 7oAf, 

in contrast to the nonrelativistic kinetic energy. 

In fact, it was a long-standing problem for the last more than 50 years in relativistic 
local field theory that the energy-weighted sum could not be expressed with the double 
commutator [201 [SB [ 22 ]. Recently, this paradox has been solved by the present authors [15]. 
The expression with commutator itself is correct, but the commutator should not be calcu¬ 
lated according to the usual way in the local field theory, where the anticommutator Eq.Q 
is employed. Eq.Q holds in the infinite momentum space. The commutator should be 
calculated first in a finite momentum space. Next, its ground-state expectation value is 
evaluated, and then we make the momentum space infinite to obtain the positive sum value. 

In other words, in the local field theory, the commutator between the time and space 
component of the nuclear four-current disappears, yielding the sum value to vanish. The 
commutator calculated in the finite momentum space, however, does not vanish, and its 
ground-state expectation value exists even in the infinite momentum space. 

Keeping this new insight in mind, let us investigate whether or not the same relationship 
as in Ea. (l35p holds in the relativistic RPA. We will first derive an expression as to the double 
commutator of the one-body operator with the total Hamiltonian, and next calculate the 
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sum of the energy-weighted strengths for exciting the RPA states, in the same way as in the 
nonrelativistic RPA [MlIIS]. 

The one-body operators, P and Q, are written as 

Here, it is not necessary for P and Q to be a function of the only coordinate. The Hamil¬ 
tonian in the a-uj model is given by 


where each term is described as 

Hq = Jd^x^p\x)ho'ip{x), (37) 

= ^ jd^xid^X2'ilj\xi)'il;\x2)Vu'il’ix2)i>ixi). (38) 

When we employ Eq.Q as in the local field theory, the expectation value of the double 
commutator with the one-body Hamiltonian P[q is expressed as 

( |[P, [Ho,Q]] \ ) = J2^<^{a\[P, [ho,Q]]\a). (39) 

a 

As mentioned before, if the operator P and Q depend on the only coordinate, the above 
double commutator vanishes. For a while, however, we keep the above form, and will be 
discussed later how the commutator should be calculated. 

The double commutator with the two-body interaction is composed of the two parts, 

( I [ -P > [ Hint , Q ] ] I )= Si + S 2 , Sk = '^ OaOa' («,«'! [ Pfc , [ Vl2 , Ql ] ] | «, a' ), 

ao.' 

where we have neglected the exchange term, and the suffixes of P and Q correspond to those 
of the two-body interaction. Employing, for example, the equation 

^Ma'(a,a'|PiVi2Qi|a,a') = ^ 0„( a | H |/3)(/3 | Rh |/3')(/3'| Q | a) 

olol' OL^0' 

= ^0„(a|PRHQ|a), (40) 

OL 

derived from the closure property in the intermediate states, the above is described in 
terms of the Hartree potential of Eq. ([5|) , 

Si = Y,B^{a\[P,[Vn,Q]]\a). (41) 

a 

The similar calculation gives 

S 2 = '^ {0^ — Ba) {0i3' - Oa') Pl3aQa'l3'Vaj3',^a'- (42) 

Finally, Eas. (l3^ . (|in) and (Il2]) lead to 

{\[P,[H,Q]] \ ) = ^Bfi{/3\[P, [h,Q]]\(3) + S2, (43) 

/3 
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where the Hartree Hamiltonian h has been given in Eq.([5|). If the first term is written, 
employing the closure property, as 

[P, [h, Q]]\P) = Y,{'^-ea)ep({l5\P\a){a\ [h,Q] |/3) 

/3 ap 

-{P\[h,Q]\a){a\P\P)) (44) 

we obtain the expression 

( I [p, [H, g]] I ) = Y,i^p-^odEo.pPpc.Qc.p 

a/3 

+ X] {^P - (^a) {6p' - 6a') PpaQa'P'Vap',pa'- (45) 

olol' PP' 


Using the notation N^p, instead of the step function, it is described as, 

( I [P, [H, Q] ] I ) = ^ (iV„^ - iV^„) E^pPp^Q^p 

a /3 

+ ^ {NaP - Npa) {Na'p' - Npia') PpaQa'P'V'aP',Pa', 

aa^/3/3' 


which is also written with the notation Ea. ()29p as 


( |[P, [P,Q]]| )= {Ppa -Pap) 

aa^/3/3' 


-^a/3,a^/3^ 

R* 4* 

'^a/ 3 ,a'/ 3 ' ^a/ 3 ,a'/ 3 ' 


Qa' 



(46) 


(47) 


In the NoN approximation, P, Q and H in the above equations should be replaced by 


P^ = 


jA^H.)APAi,(a:) , = JAi,Ha:)AQA,p(.) , 


(48) 


and 


H+ = (Po)+ + (Pint)+ , 

(Po)+= Jd^x^l;'^{x)AhoA'ij;{x), 

(Pint)+ = \ jd^Xid^X2'tpHxi)'tlj^X2)AiA2Vl2AiA2'tpix2)'lp{xi), 

with the projection operator 7l(a;) = ^ ipa{x)ipa{x). 

Ea>0 

Next, we calculate directly the energy-weighted sum, 

•Srpa = ^ a;n|(n|P|0)|^. 

n{\n = l) 

In writing the matrix elements 

{n\F\0)=^Fapd^J*, Fo,p = {a\F\f3)= f d^x ^Ux)FM^), 

ap 6 


(49) 


(50) 


the sum is expressed as 


srpa= e E 

apa'p' n{\„=l) 


(51) 
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From the relationship = F^^ and 


E /i _ \ ^ ( _,, 

n(An=l) n(An=— 1 ) 


we have 


5rpa = ^ (52) 

a^a'13' n 

Since Eos. (1261) and (152]) provide ns with 

Xn^nC^::^C^^l: = Y. + (iVa/3 " Np^) ^ 

n n ol" 

(-^q;/3 ^ “1“ (-^a/? ^f3' -^/3^Q(0^^/3^/3c 


finally we obtain the RPA sum value 

•S'rPA = - ~ Xf^a)Eap\FaY 

a /3 

^EE EafjEa'fi'iEfalS Ejja^iNa'jS' E[fjia')Vcii3\fja'■ (^5) 

a /3 a'/ 3 ' 


The first term of the above equation is nothing but the sum value in the Hartree approxi¬ 
mation, 5 h- For each approximation, it is given by. 


5'h = E Eph\Fpfi\ + < 

ph 


J2Epn\Fpn\'^ 

pn 

Y 1 Enh\Fnh\‘^ 

hn 


0 


(Full), 

(NSA, NFA), 
(NoN). 


(54) 


In the Full case, is positive, as should be so, but in NSA and NFA, the second term of 
the right hand side, coming from the N contributions, is negative owing to Fnf^ = —Fhn < 0. 
In Ref. 115]. it has been shown that the negative contributions exactly cancel the first term 
from the particle-excitations in nuclear matter, yielding Sh = 0. In writing the term in the 
Full approximation as 

pn nn hn 


n indicating the nucleon states (p -|- h), NSA and NFA neglect the first term of the right 
hand side which makes the left hand side positive. 

Now comparing Eg. (1461) with Eg. (1531) . we obtain the relationship 

5rpa = ^( |[E, [if,F]]| ). (56) 

Thus, we can express formally the energy-weighted sum in terms of the double commutator, 
as in non-relativistic RPA |16l I19j . In the NoN approximation, F and H in the above 
equations should be replaced as in Eos. ()58]) and ()59]) . 

The above double commutator, however, should not be calculated using Eq.([3|) in order 
to obtain the sum value, S'rpa, given by Eo. (l53l) . Instead of Eq.(jl|) defined in the infi¬ 
nite momentum space, we must use the anti-commutation relation in a finite momentum 
space [15]. 

{tjjaix), if^liy)} = 6abd{x - y), (57) 
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where d{x) is defined as 


d{x) = I ^ 0p = 0(Poo - IpI). (58) 

Then, we have its ground-state expectation value which does not vanish. By making Poo 
infinite in the expectation value, we can obtain the sum value. We note that Ea. ()57p is 
reduced to Eq.(jH) in the limit Poo —>• oo, because of 

d{x) (5(a;). (59) 

The double commutator in the Hartree approximation has been calculated for nuclear 
matter, according to Ea. (f57|) in Ref. [T5] . In the Eull space Na/s, the sum value of S'h is 
divergent, being proportional to P^^, in contrast to S'h = 0 in NSA and NFA. 

The energy-weighted sum for the Gamow-Teller transition strengths in the relativistic 
RPA has been explored in Ref. [23]. 


4 The continuity equation 


The current conservation in the relativistic RPA has been studied in previous papers [5l[6l 13 
|9]. One of the reasons why NFA and NSA were accepted so far is because they do not violate 
the current conservation. Generally speaking, it is important for phenomenological models 
to keep at least well-known fundamental principles. In particular. Refs. UM investigated 
the electron scattering where the continuity equation must be essential. 

There are various ways to show the current conservation in NFA and NSA. One of 
the ways is to write the Hartree polarization function in Ea. (|19p in terms of the Hartree 
propagator, Gh(1,2), 


fPH(A(I), P(2)) = Tr (aGh(1, 2)PGh( 2, 1)) , 
where Gh(1,2) is given by[5] 

1 

Gh(1,2) = — / q^(^xi,X2,uj), 

J — OO 

with 

Gh(®1,»2,w) = VvPa(a;i)^a(»2) ( -p 

\UJ — ha + l£ UJ — ha — IE 


9n 


Owing to the closure property Eq.Q, the above Hartree propagator satisfies 

- M - U{x)^Gn{x,y) = d{x - y). 

Then, Ea. dUUl) provides us with 

dxnn{l^i{x),B{y)) = -Tr (^(5(x - y) + [M + I!{x))GB.{x,y)^BGn{y,x) 

+ Tr Gn{x,y)B{5{x - y) + Gh( 2 /,x)(M -b P(x))^ 

The above equation and Eq. (II9h imply for the RPA polarization function to satisfy 

d^,n{^,{x),B{y))=d. 


(60) 

(61) 

(62) 

(63) 


= 0. (64) 


(65) 
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Thus, as far as Eg. (1631) holds, the relativistic RPA conserves the current. As mentioned 
before, Eq. ([B 5]) requires the complete set of eigenfunctions for the Hartree field. If N states 
are neglected, Ea. (l63]l does not hold as, 

- M - i:{x)^Gn{x,y) = (5(xo - yo) ^ {Pa + ha) ■ (66) 

a 

The above proof of Eq. dbhl) seems to be independent of how the Hartree ground state 
is occupied by N and N. Indeed, both NSA and NFA satisfy this equation. Physically, 
however, N states should be occupied in the ground state, as far as the negative energy 
states are required in the model framework. Hence, it may be instructive to see explicitly 
how the continuity equation is satisfied, when only a part of N-N excitations are included 
in the calculations of the excited states. 

Using Eas. dlll) to (|13p . the Hartree polarization function can be written as, 

innW),Bi2)) = (67) 

where we have used the notations: 

=Tpp{x)-i^,ipa{x)Tp^{y)B^pfi{y) 

r(“^)(xo,yo) = (NapOixo - yo) + Nf^aO{yo - xq)) 

Then, we obtain 

id^nu_{-ff,{x),B{y)) = S(xo - yo) ^ I/) “ ^/3a) 

q(/3 

+ - iEa/sB^‘^^\x,y)J t("^)(xo, yo)- (68) 

a/3 

Since Pa(x) is an eigenfunction of the Hartree Hamiltonian, the second term of the right 
hand side vanishes. Employing the expression of Eq. (| 18I) for the full space, and then the 
closure relation Eq.([3]), the above equation can be written as 

i^^^n{l^,{x),B{y)) = 5{x - y)(^ iffi{x)Biffi{x) V - X] ^c.{x)Bipa{x) ha 

^ y a 

+ '^Pj3{x)Bipp{x)nii - '^Tp^{x)B^pa{x)na \ =0. (69) 

P a ^ 

This eqnation shows the first and the second line vanishing separately. As a result, NSA and 
NFA satisfy the continuity equation, although they have the only first line. The second line 
comes from the divergent terms describing the excitations of N in the Dirac sea. Nap = {j>a + 
ha)np- This fact implies that the vacunm shonld satisfy the continuity equation by itself, 
and that the current conservation is independent of the divergence problem. Thus even if 
the continuity equation is described correctly, it is not assured that the same approximation 
is applicable to calculations of other physical quantities. 

We note finally that Eq. (|65|) provides us with the familiar form of the continuity equation 
to be expressed as the transition matrix element, 

iUnJ$a{x) ■ J pc^{x)^C^aP =0, Jpai.^) = P p{^)l^ P ■ 

a/3 
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5 The spurious state 


The present relativistic RPA equation has spurious solutions in the same way as in non- 
relativistic models [T6l [T9] . If [H ,Q]= 0, Ea. dTTl) provides us with 


PqP ) 

OLOi' PP' 




A* 

^a/3,a'/3' 



The above equation holds for any Pap^ so that we have 



-^0/3,o'/?' 
^al3,a'l3' 


A* 

^aP,a'P' 



= 0 . 


(70) 


Thus, when Qap / 0, Q is the spurious solution of Ea. (l30]l with cOn = 0. If the Eull RPA 
has the spurious state, NSA and NEA also separate it from the other solutions. As seen 
in Ea. (]29]l . A and B depend on Nap. Therefore, the spurious states can not be described 
without N-degrees of freedom not only in the Eull calculation, but also in NSA and NEA|10j. 
In the NoN approximation, [77+ , <5+] = 0 would be required for Q+ to be the spurious 
solution, but not be assured by [ 77, Q ] = 0- 

Up to this stage, it seems that the spurious state is well described in NSA and NEA. 
As discussed in ^ however, unphysical effects coming from the disregard of the divergent 
terms are hidden in their RPA equation. In order to see this fact in more detail, let us 
describe Ea. ()70p explicitly. On the one hand, it is written in the Eull case as 

PapQap T 'y ^ ip a' ^/3' T Pa'IT'P'') (y^^aP' ,Pa'Qa' P' Naa',p p'Q P '— O' (Tl) 

a'P' 


On the other hand, it becomes in NEA, 


PapQaP P y ^ Pa' ^P' ^a'^P'') {y^aP',Pa'Qa'P' Naa',pp'Qp'a'j 0, (12) 

a'P' 


and in NSA, 


PapQap T ^ y Pa' hp' T 

a'P' 



,pa'Qa' P' Naa',PP'Q P'a' 


= 0 . 


(73) 


It is seen that Eg. ([731) is the same as Eg. ([72]) . by exchanging the suffixes a' and /?' in the 
factor with fia'hp/. In writing Pahp + Pa^p = Pahp — haUp + {pa + ha)np for the Eull case, 
we can recognize in Eg. (1721) that NEA ignores the last term [pa + ha)np, so that the minus 
sign of —haUp changes unphysically the one of the interactions in NEA and NSA. 

Let us give a few comments on the above discussions. Eo. (l7Up for Eull and NSA can be 
also derived by calculating the matrix element, 

( |ata„[77,g]| ) =0 (74) 

for [77, Q] = 0. Remembering the Hartree grand state | ) to be different from each other in 
Eull and NSA, the above equation is expressed using Nap as 

NaP y ^ Qa'p' (y^a'P'i I ^p^aNcL^iCLpi \ ) Nj^ia'{ \ 0. j^CLaOj^'CLp'H \ 0, (75) 

a'P' 
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which leads to, neglecting the exchange terms, 


/ 0'Qa'P' a/3 ,a'IS'Q/3'c 


a'(S' 


= 0 . 


(76) 


Combining this with a similar equation from [77, Q^] = 0, we obtain Eq. (l70|) . 

Notice that in Ea. (|75l) is not simply replaced by the one of NFA. The step from 
Eg. ffTl]) to Eg. ([7^ requires the complete set for Na'jS'-, whereas NFA ignores the vacuum 
polarization part, {pa + ha)np, of Pahjs + Pan/s = Pah/s — hafip + {pa + ha)np in the Full 
case. Indeed, the vacuum polarization part of Eo. dTTI) does not satisfy 

^ '^j jPa' T pa'Qa'/S' ^aa'P/S'QjS'a^ — 0- (77) 

a' 


and the terms corresponding to Eo. (j72p in Eo. ()7ip also do not vanish. When Eo. (l73p for 
NSA holds, however, Eo. (l72p for NFA also does. This fact implies that neglecting the 
vacuum polarization in RPA is equivalent to assuming empty N states as in NSA with the 
complete set for In order to render Eg. (1721) valid in NFA, as well as Eg. ([73P in NSA, 
the nuclear interactions for the Full case have to be modified self-consistently together with 
Eap. Even so, as in the case of the continuity equation, the description of the spurious state 
in NFA and NSA does not provide us with a justification of their approximations. 

In the NoN approximation, [H,Q] in Eq. ifTip may be replaced by the projected one, 
[77+, Q+] = 0, in addition to the Hartree ground state. 

6 Schematic model 

It is well-known in non-relativistic RPA that a schematic model with separable interactions 
illustrates well the general character of the RPA solutions [16]. Let us explore the structure 
of NFA, NSA and the full relativistic RPA by using a similar model. 

We assume the interaction V of which matrix elements are given by 

a 


Then, Eq. (f26|) provides us with 


C^ck/3 


- NjSa 




a/3 


which lead to 


^a(3 — 


E . K 


ui — E, 


a 


1 


ajS 




iO -|- Eais 


Since we have 


■Vq ^ ^ Rab(.^')'^b i Rabija) ^ ^ 


_ ^a/3 - Npa Ja)* Ab) _ 


Q:/3 


U) 


a fa = (79) 


_ 771 J aP OijS 


the eigenvalues of the RPA equation are determined by the dispersion equation 

det KaRab(y^')^ — 9* 


(80) 
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The normalization Eg. (1311) becomes 


A = - Nf,a)\C^p\^ = - 


a/3 


In the Full RPA, Rab{(^) is given by 


ab 


dRab 

dio 


■Mb- 


flS'V) = E 

ph 

- pph 


Ja)* Ab) 
J ph ph 


oj — E, 


ph 


Aa)* Ab) ' 

J hp hp 
OJ + Epbi 


+ E 


pn 


«(a)* ^( 6 ) 
Jpn Jpn 

uj Epji 


Aa)^ .( 6 )' 
J np J np 

UJ -\- E, 


'pn 


= K» + K^A)- 


( 81 ) 


The part for particle-N excitations can be written by using Pauli blocking terms as 


with 

KU^) = E 


= «(“) - K£H, 

p(a)* p{b) ^(a)* p{b) \ 

Jnn Jnn Jnn Jn 1 


nn A 

Then, we may write 


UJ T Ejin 


= E 


A a)* Ab) 

hn ^ hn 


Aa)* Ab )' 
hh J hh 


nh 


^ Efin ^ E -E'/in 


(^) = KIA) - R^UA + RlliA- 


- pph 


hn/ 


(82) 


(83) 


(84) 


In the case of NSA and NFA, {Nap — Npa) is given by {hp — ha) in RabiA of Ea. dT^ . so 
that we have 

^NSA.NPA(^) ^ , 85 ) 

It is seen that NSA and NSA neglect R)^D(uj) in the Full expression Ea. (l84p . leaving the 
Pauli blocking term i?)j^(a;) with the minus sign. 

When the interaction has the only one component Vap^pa' = bifapfUp": have simply 


R{A = E 

a/3 


^ ~ Rap bO + Eap J 


\fpo 


A = -|V|^^. 


( 86 ) 


and the eigenvalues of the RPA excited states are obtained from the equation 


Eq. 


1 - nRioj) = 0, ^ < 0. 

doj 

becomes for the full space, and for NSA and NSA, respectively, 

iiFu„(u.) = E + E = r'-N) + «-(“). 


a ;2 - F 2 . 

ph 


- W 2 _ ^2 

pn P"' 


Rnsa,nfa{A = E 


2 Fpfe|/pfe |2 

ph 


E - R^^iA- 


nh ^bn 


(87) 

( 88 ) 

(89) 


As mentioned before, the sign of the second term of the right hand side for NSA and NFA is 
opposite to the one for the Full case. This changes the signs of k and derivative in the part 
of the N-N excitations in Eq. (l87p . The former change makes the attractive force (repulsive) 
work as a repulsive(attractive) one, and the latter change produces RPA states with negative 
excitation energies. 
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The above unphysical effects in NSA and NFA are not independent of low lying RPA 
states. When uj <C 2M*, M* being the nucleon effective mass, on the one hand, we can set 
the following equation into Ea. ()88p for the Full case. 


Ri RP®(0) 


E 


Epn 


+ R^^{0) < 0 . 


(90) 


Hence, the dispersion equation of Ea. ([871) can be written as 

«S" = 

On the other hand, in Eq. 


1 + K|RPfi(0)| 
for NSA and NFA, we may use an approximation. 


R^^{u) « < 0 , 


nh 


^hn 


(91) 


(92) 


so that we have the dispersion equation with k 


NSA,NFA 
eff 


instead of in Eq.([9l]), 


NSA,NFA 


K 

1 - K|Rhn(0)| ■ 


(93) 


In the degenerate limit where all the particle-hole energies Eph are put equal to e, the 
dispersion equation provides all the solutions but one to be trapped at the unperturbed 
energy, and the one at 

+ 2KefTe^ l/pfeP, (94) 

ph 

where Keff is given by , or 

The difference between the full space RPA, and NSA and NFA appears in the denomi¬ 
nator of KeS- The contribution from the coupling with N-N states have an opposite sign to 
each other. If the interaction is attractive (repulsive), k < 0(> 0), one should have KeS < k, 
which enhances(diminishes) the attractive(repulsive) force. This fact is realized in Eq. (l9ip 
for the full space RPA, whereas not in Eq. lfMl) for the NSA and NFA. 

In fact, these effects of the Pauli blocking terms in NFA and NSA were recognized in the 
previous numerical calculations. In the response functions to quasielastic electron scattering, 
the effects are rather negligible [6], but not in the description of the low lying states like the 
spurious state[T0] and giant monopole states [T^ [T3l ITl] . 

The RPA spurious state to be predicted at zero energy is dominated by attractive forces. 
In relativistic models, the phenomenological attractive force is rather strong. Hence, in some 
case an imaginary solution of the RPA equation was obtained, when the coupling with N- 
N states was ignored. By taking into account the coupling, the attractive interaction is 
weakened effectively, and reproduced the spurious state at zero energy [TO]. 

The same thing happened in the calculation of giant monopole states which were also 
sensitive to attractive forces. Without the coupling, the calculated excitation energy was too 
low, but the coupling played a role of the repulsive force, and explained well the experimental 
values [T2l [131E! ■ Thus, the agreement with experimental data does not mean always models 
to be physically proper for description of phenomena. 

More intuitive understanding of the difference between the full calculation, and NSA 
and NFA is shown in the Fig 1 and 2. Fig.l shows schematically the dispersion relation 
of the full RPA, Ea. (|88p . in the case of k < 0. It has a familiar structure to be found in 
literature on non-relativistic RPA|16|. The thick and thin solid curves are calculated with 
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and without the coupling with the N-N states, respectively. The solid circles denote the 
eigenvalues corresponding to A = l{dR/duj < 0), while the open circles A = —1. It is seen 
that the excitation energy of the lowest state is pushed down owing to the coupling, as 
expected. 



Figure 1; Graphical solution of the full RPA dispersion equation with Eg. (1881) . The thick 
and thin solid curves are calculated with and without the coupling with the N-N states, 
respectively. The solid circles denote the eigenvalues of the RPA excited states. The vertical 
broken lines indicate the position of the unperturbed energies. For the details, see the text. 

Fig.2 shows the dispersion relation with Ea. (|89p of NSA and NFA. On the contrary to 
Fig.l, the lowest positive energy state is pushed up with the coupling, in spite of the fact 
that the interaction is assumed to be attractive. Moreover, except for low lying positive 
energy states around Eph, the eigenvalues with A = l{dR/duj < 0) appear at the negative 
energy region around —Ehn in both NSA and NFA. 

In addition to the thick and thin curves as in Fig.l, the dotted curve is shown, which 
is obtained without the coupling and with a less number of particle-hole states. Compared 
with the thin curve, we see that the excitation energy of the lowest state is decreased with the 
increased number of the particle-hole states. Therefore, when the attractive force is enough 
for the spurious state with zero energy in NSA and NFA, its eigenvalue becomes imaginary 
in neglecting the coupling with N-N states, because of > RP*^(0) 

as in Eq. (l89p . These dependence on the number of the configurations was also observed 
numerically in Ref. [TO]. 



Figure 2: Same as Fig.l for Eq. (|89h . The dotted curve is obtained without the N-N coupling 
and with a less number of particle-hole states. For the details, see the text. 

Finally, it may be useful to re-examine in the present RPA framework the relativistic 
Landau-Migdal parameters of the a-uj model developed in Refs. ISIlTllSI. In discussing the 
response of nuclear matter at low momentum transfer g ~ 0, the interaction Ea. (l23p of the 
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a-oj model can be expressed as a separable form. It is written in Ea. dTSD as 

_ ^ / — I I /^ \ _ 

^ OLp 


/i“> = ^{a|e'‘'"r»|f() = tE2^w„r‘wi,, 


with 


' - fa = -gllml, a= -1, 

= < +fij = 9l/ml, a = 0, 

^ — fu) 1 a = 1, 2, 3. 

Here, H denotes the volume of the system which we need for rewriting the integral of Eqs.(j9|) 
and CO]) with the summation, and Wa = Ws{pcr) stands for the four-component spinor, a 
representing {s = zb, p,ct}. 


w+{pa) = 


w-{pa) = 


lEp + M* 


2Er, 


1 

ap 


Xo 


Ep + M* 


lEp + M 

2En 


* /- 


ap 


Ep + M* 


Xo 


(95) 


(96) 


1 


with Ep = -b M *2 and the 2 -component spinor, Xa- 

The matrix R in Eq.([5D]) can be divided into the longitudinal and transverse part[B|. 
Taking q = (g, 0,0), the former is the 3x3 matrix depending on a = —1,0 and 1, while the 
latter is composed of a = 2 and 3. The longitudinal part, which is required for the present 
discussions, is calculated in NFA and NS A as 


Dl = det(6ab - KaRabi^^)^ = 


sfaRlt 

-sLRl^ 


Ca + faRf faRfu. 

— f _ f pP^ 

s f PP^ s f PP^ f pPh 

where we have used the following abbreviations, 


(97) 


a = 1 - = 1 + LR, 


00 ) 


Cl = l-L[R\i + ^ 


(98) 


together with s = u/q, vp = pp/Ep^ for the relativistic Fermi velocity and A^f = 

for the relativistic density of states at the Fermi surface. Here Ep^ is defined by -b M*'^, 

using the Fermi momentum pp. The function i?P^ for isosymmetric nuclear matter is given 

by 

R'^^ = Np{l-vl)<P{x), Rt = Np^l-vl<Pix), RZ^ = Np<P{x), 


where (P{x) stands for Lindhard function with x = oj/vpq, 


<P{x) = - 1 -b - log 


1 -b X 


?7r 

-—\x\e{i-\x\ 


(99) 


1 — X 

By writing the determinant Ea. (l97p in terms of the Landau-Migdal parameters, Fq and Fi, 

as [T7] 

_ _ 1 j_ f^ pph _ foJ j^ph _|_ f‘^ „2j^ph 

El 


CaC^.Ch 


a " 

= i-(fq + 


l + i^i/3^ 


x‘^ I <^{x ), 


( 100 ) 


18 





















we obtain 


( 101 ) 


p {l-vl)Fa ^ vjF^ 

^ If phn 1 _L f phn ’ ^ if phn ’ 

J- /(t-K_i_i i + JujJ^QQ i 

with F^ = Npf^ and = lYp/o-. Thus, Pauli blocking terms yield the denominators 
depending on each meson exchange. They are calculated as, 


R 


hn 

- 1-1 


rPF d3p p2 

'o 


Fot = 0 , 


phn 

rtii 


d^p El - pV3 
Jo (2vr)3 El 


( 102 ) 


It is seen that the attractive interactions are quenched by the Pauli blocking terms in the 
same way as in Ea. (|93p . In the present model, there is no contribution to the repulsive part, 
since Fqq = 0 due to = 0. 

The parameter Fq is responsible for the nuclear incompressibility coefficient, which de¬ 
termines the restoring force of giant monopole state. Reduction of the attractive part makes 
the value of the coefficient higher [7]. 

In contrast to Fq, the parameter Fi is constrained by more fundamental requirement 
that the Femi energy Fp and momentum pp are transformed as a four-vector [25] . In the 
a-oj model, the parameter Fi comes from the longitudinal part of the w-meson exchanges 
as a relativistic effect, while the nucleon effective mass stems from the fi-meson exchange. 
They, however, are not independent of each other, as in nonrelativistic models, and should 
satisfy, according to the above requirement [11] . 


F 

^Pf 

Fp 



(103) 


for Fp = Fb -|- M, Fb being the binding energy per nucleon. As far as Ea. (ll03|) holds, 
Fi describes correctly the center of mass motion by the Lorentz boost [8], and restores also 
an abnormal enhancement of magnetic moments due to the effective mass in the Hartree 
approximation [TT] . Thus, although NFA and NSA seem to be consistent with the framework 
of the Landau-Migdal parameters, this fact does not imply that the divergence can be 
neglected. 

We note that the last term NpVp/3 in of Ea. (|98p comes from particle-hole excitations 
through the space component of the cj-meson exchange 

(104) 

Because of the last term, the continuity equation does not hold in the particle-hole space, 
-^11 7^ '*^-^00 ~ NFA and NSA, however, take into account a part of N-N states. The 

calculation of in Ea. (ll02p provides —FpUp/3, which leads to the continuity equation, 
= s^PFqq — Fqq), as discussed in 21 


7 Conclusions 

The structure of the relativistic random phase approximation(RPA) has been investigated 
in detail. The energy-weighted sum of the RPA transition strengths is expressed formally as 
the Hartree ground-state expectation value of the double commutator between the excitation 
operator and the Hamiltonian, as in non-relativistic models. In calculating the commutator, 
however, the usual anticommutation relation between the baryon fields cannot be used [15]. 
Otherwise, the sum, which should be infinite[T2|, would vanish. 
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The main difference of the relativistic RPA from the nonrelativistic one stems from 
antinucleon(N) degrees of freedom, but they cause the divergence problems. The two kinds 
of approximations were proposed by previous authors in order to avoid the problems without 
the renormalization. The one[5l |6] is the no free term approximation(NFA) which simply 
neglects the divergent terms in RPA response function, and the other [TO] is the no-sea 
approximation (NS A) where N states are assumed to be empty. Actually, both approxima¬ 
tions are equivalent to each other. They were employed widely and shown to work well for 
reproducing experimental data in a phenomenological way[3l E] El [TOl IlSl [TO] [14]. 

The present paper has shown that NFA and NSA have the serious problems. The RPA 
dispersion equation yields the RPA states with negative excitation energy, in addition to 
the low lying positive energy states. This fact implies that the RPA ground state is not the 
lowest one. Owing to those negative excitation-energy states, the energy-weighted sum of 
the transition strengths vanishes. These results are not avoidable for NFA and NSA which 
satisfy the RPA relation of the energy-weighted strengths, since the relativistic sum rule 
value stems from the excitations of Dirac sea[T5]. Moreover, since the only limited space of 
nucleon-antinucleon states is included in NFA and NSA, attractive(repulsive) forces work as 
repulsive(attractive) ones between their couplings. This fact affects also the couplings of the 
particle-hole states with nucleon-N states. Unfortunately, these unphysical couplings played 
an important role in explaining the spurious state[TO] and the giant monopole states [TOl [TOl 
E] in the previous numerical calculations in NSA. These results have been shown clearly 
by using a schematic model. 

R has been shown that there is no problem for NSA and NFA to describe the continuity 
equation, since it is independent of the divergence. 

Thus, N degrees of freedom which provide the divergence are not ignored. As far as a 
part of the N space is necessary, the rest of the space also should be taken into account in 
a proper way, even in phenomenological models. Indeed, it was shown in Refs. laiTllSSllST] 
that the renormalization of the divergence plays an important role in discussions of some 
physical quantities. Those roles are state-dependent, and could not be incorporated into 
phenomenological interactions or their coupling constants. Moreover, if the divergence of the 
linearly energy-weighted sum is understood, we can make clear the meaning of the analyses 
as to the distribution of transition strengths with the energy moments]!]. 
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